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£f) ■ Abstract 

We study the behavior of the antiferromagnetic RP 2 model in d — 3. 
The vacuum structure is analyzed in the critical and low temperature 
\ regions, paying special attention to the spontaneous symmetry breaking 

■ pattern. Near the critical point we observe a full breakdown of the 0(3) 

symmetry of the action. Several methods for computing critical expo- 
nents are compared. We conclude that the most solid determination is 
obtained using a measure of the correlation length. Corrections-to-scaling 
|' are parameterized, yielding a very accurate determination of the critical 

D , coupling and a 5% error measure of the related exponent. This is used to 

i-C ' estimate the systematic errors due to finite-size effects. 
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1 Introduction. 



Antiferromagnetic, (AF) interactions give rise to very interesting properties in 
statistical systems. Specifically, the symmetries of the broken AF phase can 
be very different from their ferromagnetic counterparts. Ground states with 
frustration or disorder are common in these systems. 

oThe interest of AF models extends from the study of some physical prob- 
lems where the antiferromagnetism appears naturally to the more theoretical 
aspects of new Spontaneous Symmetry Breaking (SSB) patterns, or Universality 
Classes. Moreover, attention has been recently paid to AF models in four di- 
mensions [0, ||. There, one could even hope to gain some insight on the puzzling 
problem of the non perturbative formulation of interacting field theories. 

In this paper we will be concerned with the RP 2 model in three dimensions. 
Its ferromagnetic sector has been extensively used to study liquid crystals ||. 
The AF one has been also studied |I| ||, |) as it presents an unfrustrated ground 
state with nonzero disorder. It seems to belong to a class of models with a 
SSB pattern of type SO(3) x SO(2)/SO(2) In condensed matter physics 
it is not rare to find systems with an equivalent SSB pattern, like super-fluid 
He , helical || and canted || spin systems and some frustrated quantum AF 



Heisenberg models |10 , the latter being specially interesting because of their 
possible relation with High Temperature Superconductivity JHJ . 

A radical reading of the Universality Hypothesis, suggests that the critical 
properties of a model are completely fixed by the space dimension and by its 
low and high temperature phase symmetry groups, H and G respectively. Even 
more, it is also generally assumed that models with locally isomorphic mani- 
folds G/H have identical critical behaviors, regardless of the global properties. 
This scenario has been theoretically challenged in ref. JO, where it has been 
suggested that the massive modes, not fixed by the local properties of G/H, 
could change the numerical values of the critical exponents. 

The above SSB pattern has been studied in perturbation theory Jl^j , where 
the main conclusions reached are that the only possible critical points in three 
dimensions are second-order, either with the Mean-Field (MF) exponents or 
with those of the 0(4) sigma model, or first-order. 

AF interactions frequently produce first-order transitions jli], |j, although 
with a large correlation length (weak-first-order transitions) that are charac- 
terized by the (apparent) critical exponents: v = 0.5, 7/f = 2, a/v = 1 JTsf . 
These exponents are found in ref. (l6| for the Stiefel manifold V^. However, 
other Monte Carlo (MC) simulations for helical and canted spin models ^ [TtJ . 
have yielded exponents not far but incompatible with those of a weak-first-order 
transition. 

On the other hand, the RP 2 antiferromagnetic model has a critical point 
with clearly different from weak-first-order exponents {y « 0.78, ajv ss —0.44) [|| 
and very close to those reported for 0(4) jl8|, although hardly compatible with 
them. 

The actual relevance of the model in experimental situations is difficult to 
know. It is an uneasy task to compare the critical exponents of this model 



2 



with experimental data. In fact, the results for spin systems are not very clear, 
even regarding the order of the phase transitions encountered (see ref. || for a 
review). For the super-fluid 3 He, the possibility of facing our critical exponents 
with experimental measures is remote, the critical region being so small that it 
seems experimentally inaccessible [ fL9[ . 

In this paper we present a more detailed exposition of the studies of ref. Q , 
completed with new analysis of the MC data. Our statistics have also increased 
by an amount of about a 50%. The high statistics has allowed a thorough 
study of the autocorrelation times, both exponential and integrated, for sev- 
eral observables. We are able to parameterize the corrections-to-scaling. As a 
consequence, we get a high precision determination of the critical point, and 
we measure the first correction-to-scaling exponent with an error of 5%. We 
use this exponent to extrapolate the critical exponents to their infinite-volume 
values, and get an estimation of the systematic errors. In addition, we carry out 
a study of the ordered phase both numerically and analytically. We find a low 
temperature phase with a global 0(2) symmetry and, from finite-size scaling 
analysis, evidence of a breaking of this symmetry close to the transition. The 
MF approximation is considered at the critical point which already shows the 
extreme complexity of the symmetry-breaking pattern in that region. We also 
analyze the continuum limit of the model obtaining the SO (3) x SO(2)/SO(2) 
non linear sigma model. 

We start in section ^ by defining the model and the observables, and discuss 
some theoretical aspects and approximations. The MC algorithms used and the 
dynamics of the simulation are described in section In section [I| we present a 
determination of the critical exponents using finite-size scaling techniques. The 
reliability of a parameterization of the corrections-to-scaling is considered in 
section [|. Finally sections || and [?] are devoted to discuss the vacuum structure 
of the broken phase. The conclusions are in section H. 



2 The Model. 

We will consider a system of spins {vi} taking values in the sphere S 2 C R 3 and 
placed in the nodes of a cubic lattice. The interaction is defined by the action 

<i,3> 

where the sum is extended over all pairs of nearest neighbor sites. The partition 
function is constructed as 




where the usual 0(3)-invariant measure is used. 

The action and the measure are invariant under a local Z2 transformation 
(vi — > — Vi). As this is a local symmetry, Elitzur's Theorem guarantees its 
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preservation even after a SSB. This means that the natural variable for the 
model is in an equivalent class of the quotient group RP 2 = S 2 /Z2, the two- 
dimensional real projective space. The action (|l|) can be rewritten in terms of 
the RP 2 variables, which in turn can be expressed as the tensorial product of 
the original spins Tj = Vi <£> Vi, 

S = -0^2tvr i T j . (3) 

<ij> 

As the measure for the t is much more complex than for the vectors, it is 
generally more convenient to work directly with the latter. 

2.1 Observables and measures. 

The observables to be measure in a Monte Carlo simulation on a finite lattice 
must be scalars under the 0(3) global symmetry group and obviously also under 
the local Z2 group. Let us construct the traceless tensorial field Tj, whose 
components are 

Tf =v?vf- V. (4) 
The Fourier transform in a L x L x L lattice can be defined as follows: 

r£L 3 

where we denote by r the position vector of the site, and by p the momentum 
(p = (2im x /L, 2im y /L, 2irn z /L), with rij = 0, ...,L — 1). 

In principle, we are only interested in the magnetizations associated with a 
fixed momentum. In addition to the (normalized) zero momentum magnetiza- 
tion 

M = y Zl T (W) = y T (0,0,0), (6) 

x,y,z 

it is interesting to study the staggered magnetization related with states of 
period 2 under translations, 

M S = ^E (-l) X+ * +Z T (Kl3/ ,,) = iT (w) , (7) 

x,y,z 

that is the (normalized) difference between the magnetization of the even sites, 
defined as those where x + y + z is even, minus that of the odd ones, defined 
analogously. 

Other finite momenta (nonzero in the L — > 00 limit) could be necessary if the 
vacuum would present higher period symmetries. However we have not found 
any reason, neither analytical nor numerical, to expect this. 

We compute the magnetization and the susceptibility respectively as 

M = (VtrM 2 ^ , (8) 
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X = V(trM 2 ). 
It is also useful to consider the quantity 

_ ((trM 2 ) 2 ) 
K (trM 2 ) 2 ' 



(9) 



(10) 



which is directly related with the Binder parameter [ f20[ . 

We use the second momentum correlation length because it is easier to mea- 
sure than the exponential (physical) one, but it is expected to have the same 
scaling behavior at the critical point [Ell E2| , 



X/F-I 
4sin 2 (7r/i) 



1/2 



where F is defined, using the shorthand notation |T| 2 = TTT, as 



-(ttAO) 



+ 



-(0,^,0) 



-(o,o,¥) 



(11) 



(12) 



Analogously we define the associated quantities for the staggered observ- 
ables. 

We also measure the nearest neighbor energy 



<i,3> 



and the next-to-nearest neighbor energy 



(13) 



(14) 



<l,3> 



which is useful to study the 0(3) broken phase. 

The measures of E\ have been also used to calculate the (3 derivatives 
through its connected correlation with every observable. 



2.2 The ground state. 

When j3 ^> the ground state corresponds to a configuration in which all the 
spins are aligned in a given direction. This state is thus 0(2) invariant. The 
thermal fluctuations do not destroy this ordering until /? ~ 1.3, where a first- 
order phase transition occurs. This system has been used for modeling liquid 
crystals 0. 

For negative (3, the ground state corresponds to a highly degenerated non- 
frustrated state where all spins are orthogonal to their neighbors. 

The zero energy configurations that dominate the measure are those in which 
the even (odd) sublattice is aligned, and the spins in the odd (even) one lie 
randomly in the orthogonal plane. It is easy to check that other subsets, such 
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as, for example, those with spatial periodicity 3 Q or without any periodicity, 
have a relatively null contribution to the measure. 

When —j3 is large but finite, the fluctuations effects must be taken into 
account. Let us consider a spin of the aligned sublattice surrounded by its 
neighboring spins lying in the orthogonal plane. Should they prefer to orientate 
in a given direction of the plane, that tendency would permit a larger fluctu- 
ation of the spin of the aligned lattice. As a consequence, there is an induced 
interaction between second-neighbor spins in the plane-aligned sublattice. One 
can think of a RP 1 model (that is equivalent to an 0(2) model) in a face cen- 
tered cubic lattice with a nonzero coupling. The important question is if the 
induced coupling is large enough to produce an 0(2) phase transition. The 
Monte Carlo results show that, close to the antiferromagnetic phase transition, 
the 0(2) symmetry is also broken, but it is restored for lower values of (3. We 
address to sections || and [?] for a detailed discussion on this subject. 

2.3 The Mean-Field Approximation. 

A standard tool in statistical mechanics is to study the system, close to the 
critical point, by means of the MF approximation. This analysis is known 
to improve as the system dimensionality grows. It is also useful to know the 
predictions of a MF calculation since any deviation may be interpreted as a sign 
of non-triviality. 

We shall use the formulation of the MF approximation as a variational prin- 
ciple (see for instance p3|). It can be stated as follows: let us consider a con- 
figuration space, and the Boltzmann measure D/u = (ZfT))" 1 exp(7i(/i)/T)<i/i, 
where T is the temperature. For any other Hamiltonian 7i p (/i), the following 
inequality holds, 

T <T P + (H-H p ) P = g p , (15) 

where by T we mean — TlogZ(T), by (0) p we refer to the mean value of 
the operator O, calculated with the Boltzmann measure of the Hamiltonian 
Hp(n), and T p is the free energy correspondent to H p (/Lt). 
To study the action (|l|) we have used the Hamiltonian 

i even i odd 

where Q c ,o are tensorial mean fields. The usual procedure consists on choosing 
the parameters of those fields in order to minimize Q p . For the sake of simplicity, 
we restrict ourselves to the case of commuting Q . o - It is easy to show that in 
this way Q p depends on four parameters in addition to the coupling. The mean 
field equations are very complex, but they can be simplified in the neighborhood 
of a continuous transition. 

We report here the qualitative results addressing the reader to appendix 
for details. Regarding the critical exponents, we find the expected results 
a = and (3 S = 1/2, where the latter is the exponent related to the staggered 
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magnetization. However, the exponent associated with the non-staggered mag- 
netization turns out to be (3 — 1. Other interesting result regards the structure 
of the ordered vacuum. Performing a Ginzburg-Landau expansion of Q p on M s 
we observe that, up to fourth order, an 0(2) symmetric vacuum has the same 
free energy as a broken one. The sixth order term favors the unbroken vac- 
uum. However, this quasi degeneracy suggests that the problem of the 0(2) 
breakdown is beyond the reach of the MF approximation. 

2.4 The continuum limit of the 0(2) broken model. 

In reference p3| the continuum limit of systems with a SSB pattern of type 
S0(3) x SO(2)/SO(2) was considered with the conclusions regarding the order 
of the transition mentioned in the introduction. 

We have found that as (3 — > — oo (see section [?] for details) the vacuum is 
0(2) symmetric, and it cannot be described with the above mentioned SSB 
pattern. However, in the critical region, the remaining 0(2) symmetry seems 
to be broken, and then, we should consider a pattern of type 0(3)/{l}, which 
is locally isomorphic to S0(3) x SO(2)/SO(2). Although in the critical region 
one cannot assume infinitesimal fluctuations, it is possible to relate our model 
with one in which we add an explicit (large enough) ferromagnetic interaction 
between second-neighbors. That model presents a vacuum with the same sym- 
metries as the one with the action ([!]) in the broken critical region. It can be 
shown (see appendix B for details) that, in the continuum limit, the action of 
this model can be written as 

S= /d 3 a^[P(i? T (V?) 2 ], (17) 
J n 

where R(x) is an S0(3)-valued field and P is a diagonal matrix of couplings of 
type {gi, 51,-52} • 

The action ( |l7j ) is in the class studied in ref. 

2.5 The low temperature effective model. 

For —f3 very large, it is useful to consider a limiting model as follows. Let us 
suppose that the fully aligned sublattice, say in the z axis direction, is the even 
one. Writing the even and odd spins as 

v c = (v x ,v v ,^/l - (v*) 2 - (vv) 2 ), (18) 
v° = (v 7 ! - (v z ) 2 cos(p, y/l - (v z ) 2 sm<p,v z ), (19) 

respectively, we obtain the approximate action (a = x, y, z) 

S = -P ]T « cos^- + v\ sin^j + v]f + 0{{v a ) 4 ). (20) 

<.i,j> 
i even 
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Rescaling the spins as (3{v a ) 2 — > (v a ) 2 we conclude that in the large — (3 limit 
the tp's correlation functions become /3-independent. 

The action (|2^) is easy to simulate, although one should be aware that the 
alignment direction in a finite lattice rotates (Goldstone modes). To avoid this 
problem, a global rotation after every Metropolis sweep should be performed. 
We are particularly interested in the correlation of the second neighbor angles: 
(cos 2 ((fi—(fij)) — E2. In section]?] we check numerically that this model smoothly 
joins with RP 2 at large —f3. 



3 MC simulations. 

We have used a standard Metropolis algorithm for the update. The fluctuation 
of the spins in the critical region is very large and one can choose as change 
proposal a value almost independent of the initial state. In fact, an uncorrelated 
spin is accepted with a 30% of probability. We have used an almost independent 
proposal but with 3 hits to obtain about a 70% of mean acceptance. 

We have also checked the efficiency of cluster algorithms (24J as the model is 
suitable for applying the embedding procedure p5| . Unfortunately, as in most 
antifcrromagnetic systems, the usual cluster methods do not reduce the expo- 
nent z for the autocorrelation time (AT). We have studied both the Swendsen- 
Wang's method and the Wolff's single cluster version. The system always 
presents a cluster with around a 65% of the spins, and the remaining 35% 
forms very small clusters (only about a 1% contain more than 10 spins). The 
results regarding the efficiency are slightly worse than those from the Metropolis 
method. 

The simulations have been done distributed over several computers based on 
ALPHA, SPARC and Pentium processors. The total computer time employed 
has been the equivalent of 18 months of ALPHA AXP3000. We measure every 
10 sweeps and store individual measures to extrapolate in a neighborhood of 

The 



the simulation coupling by using the spectral density method |26j, [27 
number of performed sweeps in the critical region is displayed in tablejl . Some 
other shorter runs have been done at the peaks of the connected staggered 
susceptibility, and in the broken phase in order to study the vacuum. 

Let us first consider the statistical quality of our data. Following the notation 
of reference , we compute the unnormalized autocorrelation function 

N-t 

C {t) = J ^— t Y J O l O l+t -^ (21) 

i=l 

where Oi are successive measures for the operator O, fio is the mean value of O 
and N is the number of measures. We define also the normalized autocorrelation 
function as 
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The statistical error in the measure of O is proportional to \/2r^ t /N where 
the integrated AT Tq* can be obtained from 

T$ t (t) = ± + Y / Po(t') (23) 

i'=l 

for large enough t. In practice this value is selected self-consistently, usually as 

t = 67§*(t). 

In table |l| we summarize the results of the simulations at j3 — —2.41. In 
the L = 24 lattice, however, two thirds of our data come from a simulation at 
= —2.4. We observe that the integrated AT for the staggered susceptibility 
(with t = 6T"'*) grows quadratically with the lattice size, as expected [z ex- 
ponent equal to 2). The AT for the other operators in the table are smaller. 
We should point out (see fig. |l|) that we need a larger window to obtain a sta- 
ble value for those cases. Specifically we have used t — lOr™* and t = 20?^* 
respectively. 

As the number of measures that we have performed is very large compared 
with the AT (10 5 r£f for L = 32 and 10 4 t^ for L = 48), we can discard some 
hundreds of r™* for a safe thermalization, and make bins of consecutive data of 
that size to assure that they are uncorrelated. 

However, it is interesting to look at the exponential AT to know if there is 
some information that remains at time scales larger than the integrated AT. 
As the exponential AT measures the interval that the system remains out of 
equilibrium along a MC simulation, only when the number of MC iterations is 
much greater than the exponential AT, we can have confidence on the obtained 
measures. We do not need to assume, as is usually done, that the exponential 
AT is of the same order of the integrated one, as our statistics is good enough 
for a reliable estimation of the latter. 

The exponential AT for the operator O can be measured as the t — ► oo limit 



Table 1: Number of Monte Carlo sweeps performed for different lattice sizes. 
Measures have been taken every 10 sweeps. The integrated AT (in sweeps) for 
both magnetizations and for the energy are also displayed. We have discarded 
in each case about 200t Xb iterations for thermalization. 



L 


MC sweeps(xlO ti ) 




T x 




6 


6.71 


7.37(3) 


5.81(2) 


6.02(2) 


8 


17.07 


11.41(4) 


7.41(3) 


7.30(3) 


12 


6.51 


24.9(2) 


13.23(12) 


10.68(11) 


1G 


22.14 


44.5(3) 


22.27(14) 


14.41(8) 


24 


8.77 


107(3) 


51.3(10) 


24.7(6) 


32 


28.51 


175(6) 


90.0(8) 


34.8(3) 


48 


3.93 


410(20) 


222(10) 


60(3) 



9 



20 



15 



10 



5 







1 1 1 1 1 1 1 

~_ ° 


i i | i i i i | i i i 


1 1 1 


1 1 


.:/ . . ■ ■ ■ 








■ X 

- ' / E, 


X 












i i i i i i i 


1 1 1 1 1 1 1 1 1 1 1 


E l 

, 1 , 


1 1 



20 40 60 80 

t 

Figure 1: Autocorrelation time (in measures) estimators for Xs, X an d E\ in 
a L = 32 lattice at j3 = —2.41. The dashed line correspond to the integrated 
time, the solid one to \ and the symbols to (circles, squares and 

diamonds for Xs, X an d E\, respectively). 



We should take the supreme of these quantities for all operators. 

Usually it is not easy to find a clear t — > oo limit, since for finite t there is a 
systematic error due to the faster modes contributions. In that case, it is useful 
to study the (very noisy) estimator 

^o P{2 Ht)=lo g -^-. (25) 
po[t + 1) 

In fig. [l| we show t^*, t™*, and Tg* in L = 32 in the critical region. In all 
cases a plateau is clear. The estimators j^ xp ( 1,2 ^ for all the three quantities are 
also plotted. We observe that for Xs the three estimators are almost equal. For 
the other two cases, there is no a visible plateau for Tq Xp( ' 1 ' ) , but the estimator 
r cxp(2) seems £ Q stabilize at a value near that corresponding to Xs- We remark 
that just with the data from E\ or x we would conclude that the AT for this 
lattice size is about 200 sweeps, that is about 6 times the integrated AT for E\. 
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Figure 2: Statistical errors of \ s (upper curve) and E\ (lower curve) as a function 
of the bin size (in measures) for L = 32. The vertical dashed line corresponds 
to the bin size used in the analysis. 

A further check of statistical independence can be done by studying the 
errors computed with bins of data of increasing size. This could allow us to 
observe AT at scales much greater than the previous. 

In fig. H we show the statistical errors for the energy and the staggered 
susceptibility in L = 32 as a function of the length of the bins, using all the 
statistics. When their size is not big enough, the data in different bins are 
correlated and so, we see a growing error. When the size of the bins is enough 
to consider them independent, there is a plateau, that indicates the correct 
value for the statistical error. We can estimate the integrated AT from the 
quotient between the v alue of the error at large Lbin and at Lbin = 1, since 
c(ibin)/o'(l) —> V2r int . Using this procedure we obtain rj?* = 172 sweeps and 
T Ei = 34 sweeps, in good agreement with the values shown in table [|. 

From these figures we also observe that our choice of the number of bins, 50, 
is completely safe, allowing an accurate estimation of the statistical errors. We 
have chosen a large number of bins to ensure a 10% precision in the statistical 
error determination. 

From our results we find very unlikely the existence of a much greater AT 
for any other important observables. 
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4 Measures of Critical Exponents. 



Recalling that this model has two different order parameters (q^) and, in prin- 
ciple, two different channels for the correlation lengths, it is compulsory to 
study separately both critical behaviors. Therefore, in addition to the usual 
exponents (v,/3,j and if), we consider the analogous for the staggered channel 
(z^ s ,/3 s ,7s and ?y s ). Although the most economic scenario requires v = v s , in order 
to consistently define a single continuous limit, we shall check this assumption. 

4.1 Finite-size scaling. 

Our measures of the critical exponents are based on the FSS ansatz. Let 
(0(L,P)} be the mean value of an operator O measured on a size L lattice, 
at a coupling value /3, and let £(L,/3) be a reasonable estimator for the corre- 
lation length on a finite lattice, such as the one shown in eq. (ffTf). Then, if 
0(oo,/3) ~ |/3 — P c \~ x ° , from the FSS ansatz one readily obtains P9[ 

(0(L, (3)) = L x °l v [F (t(L, ff)/L) + L-»G (Z(L, P)/L)] +..., (26) 

where u) is the universal exponent associated to the leading corrections- 
to-scaling. The dots stand for higher order scaling corrections and £(oo)~" 
terms, negligible in the critical region. The above expression is straightforwardly 
generalized to functions of mean values. 

We remark that to obtain eq. j26| ) the condition on the definition of the 
correlation length is that £(L, (3)/L should be a smooth monotonous function of 
L/£(oo, (3) and that x% = v. That condition is verified by the quantity defined in 
eq. ([ll]) but also by kL (see eq. (|To|)) and their staggered counterparts. However, 
the scaling function Go is not invariant under a change of the correlation length 
definition, and a proper choice can largely reduce it. 

Let us emphasize that the fulfillment of the hyperscaling relations is an a 
priory condition of the ansatz. In fact, if a relation like ( p6[ ) holds, for example, 
for the magnetization and its square, one readily obtains 7 + 2/3 = dv. However, 
we will check the observance of the hyperscaling relations as a test of consistency. 

4.2 Measuring at the maxima. 

Usually, the application of FSS consists on measuring the mean value of an ob- 
servable for different lattice sizes in the apparent critical point, where L/£(L, (3) 
is assumed to be constant, so the exponent of L is easily extracted, after ensur- 
ing that the corrections-to-scaling are small. Typically one can define two such 
points, one is the infinite volume critical point (see section ^|), the other is the 
peak of some quantity like the connected susceptibility defined as 

X c = X - VM 2 . (27) 
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Figure 3: Connected staggered susceptibility as a function of lattice sizes (from 
top to bottom) L — 24, 16, 12, 8. The left limit of the (3 axis corresponds to (3 C . 

The advantage of measuring at a peak is that the apparent critical point is 
self defined, and also that the statistical error is much smaller at a peak (zero 
(3 derivative) than when the slope is large. 

In fig. H we plot the connected staggered susceptibility for several lattice 
sizes. The most significant property is the large shift of the apparent critical 
point, that makes necessary to perform simulations at values of the couplings 
dependent on the lattice size, as (3 C cannot be reached with the spectral density 
method. The MC data used in fig. || correspond to short runs performed at the 
coupling values signed by filled symbols. 

The shift of the (3 value corresponding to the maximum of the staggered sus- 
ceptibility, Afic , can be used to measure the thermal exponent v and the critical 
coupling in the thermodynamic limit. The expected behavior is A/J^ oc l^' v . 
Performing a three parameters fit, we obtain for L — 8, 12, 16, 24: v s — 0.83(11), 
/3 C = -2.410(22), with x 2 / d -°- f = 0.05/1. 

We can also compute r] s from the value of the susceptibility at the maxi- 
mum. The fit gives r) s = 0.035(39), with x 2 /d.o.f. = 0.98/1. The value for r) s 
is compatible with those obtained below, but the errors are 30-times greater 
(the statistics used is 200-times smaller). This method has three important 
drawbacks compared with the methods discussed below: a) it is less precise; b) 
it needs simulations performed far from the critical point at a coupling value 
that depends on the observable to study; and c) large corrections-to-scaling are 
expected due to the large distance from the critical point. 

Other quantities, as (3 derivatives of magnetizations or correlation lengths 
also present peaks, but they are usually very broad and noisy what implies a low 
quality of the obtained critical exponents. The specific heat does not present a 
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critical divergence as the a exponent is negative (see subsection |5.2|). 



4.3 The quotient method. 

Fortunately, one can use another method that does not require a previous knowl- 
edge of the critical point, and which is directly based on eq. (|26|). Let O be an 
observable whose mean value is measured at the same coupling in two lattices 
of sizes L\ and L2 respectively. Using (Eq) we can write for their quotient 



where s = L2/ L\. In the case of O = £, the exponent x^ is just v. Choosing 
the coupling in such a way that Q{ = s we can directly compute the exponent 
for any operator as 

Qo\ Qt= . = a? o/v +0{L- u ). (29) 

As the exponents are obtained from just two independent simulations, it is 
possible to do a very clear statistical analysis. The use of the spectral density 
method avoids the necessity of an a priori knowledge of the coupling where the 
quotient of correlation lengths is the desired. 

Let us show with an example how does it work. In fig. ^ (upper part) we 
plot the quotient between the staggered correlation lengths for L 2 /Li = 2 as 
a function of the correlation length in the smaller lattice, in lattice size units. 
For L\ large we should obtain a single curve. The deviations correspond to 
corrections-to-scaling. If we plot the quotient for the staggered susceptibility, 
we observe similar corrections-to-scaling (see fig. [| middle part). Plotting Q Xa 
as a function of Q^ 3 the corrections-to-scaling are strongly reduced even for 
lattice sizes as small as L\ — 6 (see fig. [|, lower part). 

With this technique, we have studied several operators as the susceptibility 
(x x = 7), the magnetization (xm = —0) as well as their staggered counterparts. 
To compute the v exponent, we use (3 derivatives of observables. For example, 
x d£,/d0 = v + !• I n the next two subsections we shall separately study the 
determination of v and rj type exponents. 



4.3.1 The exponent v. 

In the first two columns of table ^[ we show the values obtained for the v 
exponent measuring, at Q^ s = 2, the j3 derivatives of £ s and £. The errors 
have been computed using the jack-knife method. The corrections-to-scaling 
are smaller than the statistical errors, even for lattice sizes as small as L = 6. 

We have repeated the analysis, using the non-staggered correlation length as 
the variable. In columns 3 and 4 we write down the obtained results. These are 
fully compatible with the previous, but with a slightly greater statistical error. 

Finally (columns 5 and 6 of table ^|) we have alternatively computed the v 
exponent using the logarithmic derivative of the magnetizations, which scales as 
L x l v . The results have a smaller statistical error, but they suffer from important 
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0.48 0.50 0.52 0.54 




Figure 4: Quotient of staggered correlation lengths for L^/Li = 2 as a function 
of £ s (£i)/£i (upper side). In the middle we plot the same quantity for the 
staggered susceptibility. In both cases the higher line corresponds to L\ = 6 
with lattice sizes increasing with decreasing heights. Finally (lower side) Q Xb as 
a function of Q^ s for L2/L1 = 2 is plotted. 

corrections-to-scaling, and it is necessary to compute the L — > 00 extrapolation. 
This introduces uncertainties reducing the quality of these numbers as compared 
with the previous. Other quantities such as /3-derivatives of Binder's Cumulants 
suffer from statistical errors about three times larger, and therefore are quite 
useless for comparison. 

We remark that all the determinations of the v exponents (regardless the 
observables used being staggered or not) are compatible within errors (at the 
1% level). One relevant question is whether summing data from different ob- 
servables it is possible to reduce the statistical error, or even if the differences 
between columns are significant. As all the used quotients are flat functions 
of Q^, the main contribution to the error is due to the observable (not to the 
correlation length). We have studied the mean and the difference of the first 
two columns. The reductions of the error bars in the first case are negligible 
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(less than a 10%) because the data are strongly correlated. In the case of the 
difference, the correlation makes the error slightly smaller (about a 10%) than 
each term. The results do not show any significant deviation from zero. 

The results for pairs with s — 3/2 or s — 4/3, that can be formed with 
the simulated lattice sizes, are compatible. As our better statistics are in the 
L = 16, 32 lattices, we only report the results for s = 2. 

4.3.2 i] exponents. 

The r\ exponent can be obtained from the scaling of the magnetization or of the 
susceptibility, using the scaling relations: 2/3/v = d — 2 + rj and j/v = 2 — -q. 

The results for the ?7 S exponent are given in table |^. They have been obtained 
studying, in the staggered channel, the susceptibility, the magnetization and the 
maximum eigenvalue of the magnetization tensor (that scale with 7 S , /3 S , and (3 S 
respectively). Technically, the computation is somewhat different than in the 
case of the v exponent, since the quotients change now very fast with Q^ s . One 
has to consider carefully the statistical correlation of the data. Thus we are able 
to measure ^a/v with an accuracy of about a 0.1%, from which we can compute 
?y s = 2 — 7 S /^ with a 5% error. 

For the non-staggered channel the results are reported in table [|. The 
columns 2, 3 and 4 correspond to the use of Q^ B as the variable. We observe 
important corrections-to-scaling. The use of does not improve the results. 
Using Q k l, we find a quite strong reduction of the corrections-to-scaling. 

5 Corrections to scaling. 

For the determination of the corrections-to-scaling exponent, lo, we have used 
the property shared by different observables (Binder type parameters, k, and 
correlation lengths divided by the lattice size) that asymptotically should cross 
at the critical point. It can be shown [^0| that the scaling functions correspond- 
ing to lattice sizes Li, L2 cross at a (3 value whose shift from the critical value 



Table 2: v exponent obtained from a FSS analysis using data from lattices of 
sizes L and 2L. In the first row we show the correlation length that we ask for 
duplicate, and in the second one the operator whose quotient we measure. 



L 


d£s/d/3 


dC/d/3 


d£ s /d/3 


d£/d/3 


dlogAf s /d/3 dlogM/d/3 


6 


.786(6) 


.790(6) 


.774(14) 


.759(16) 


.689(4) .670(4) 


8 


.785(4) 


.781(4) 


.771(7) 


.774(7) 


.717(3) .706(2) 


12 


.789(8) 


.782(9) 


.781(18) 


.781(15) 


.739(9) .731(8) 


16 


.786(6) 


.780(5) 


.782(8) 


.781(6) 


.755(4) .751(4) 


24 


.77(2) 


.77(2) 


.77(3) 


.77(2) 


.756(16) .758(13) 
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behaves as 

A(3 Ll ' L2 oc 1 7 S U Li U ~ h , (30) 
S" — 1 

where s = L%jL\. In fig. ^| we plot the values of the coupling where several 
scaling functions cross, as a function of L 2 . We use along this subsection the 
value v = 0.783. According to ( |30| ) if we fix L\, the asymptotic behavior should 
be linear, what is well satisfied in the plot. 

We have performed a global fit of the points from each observable to a 
functional dependence of the type of eq. ( (30| ) with L\ — 6. The good precision 
obtained for the v exponent makes negligible the uncertainty that comes from 
this quantity. In fig. ^| (upper part) we plot the values of the extrapolated (3 C from 
independent fits for every observable after fixing the value of the corrections-to- 
scaling exponent lu. 

Studying the variation of the % 2 function (using the full covariance matrix) 
we obtain 

k s : /3 C = -2.4097(19), lu = 1.05(45), 

f/L : f3 c = -2.4096(23), lu = 0.75(20), (31) 

\ a /L : p c = -2.4082(13), lu g (0.6, 1.9). 

The corrections-to-scaling in the k case are so small that it is not possible to 
measure the parameters. Repeating the fits with L\ = 8 we obtain very similar 
results. 

From these fits we conclude that all the extrapolated values coincide up to 
A/3 = 0.002. Notice that this implies that both correlation lengths diverge at 
values of (3 which are at most 0.002 apart. We can go one step further under the 
hypothesis that the extrapolated values must be the same, which imply that the 
infinite volume propagator has two separated poles in momentum space, one at 
the origin, the other at (n/a, it /a, ir/a). A x 2 minimization of all the (strongly 
correlated) data gives (see fig. [| middle and lower parts) 

P c = -2.4088(3), lu = 0.87(3), x 2 / d -°-f- = 19.3/14. (32) 
We have repeated the fits discarding the L\ — 6 data, and measuring the crossing 



Table 3: r/ s exponent obtained from the FSS of the staggered susceptibility, 
magnetization and maximum eigenvalue of the magnetization, using £ s as the 
variable. In all cases, the pairs are of type (L, 2L). 



L Xs Ms K 



6 


0.0431(10) 


0.0474(9) 


0.0480(12) 


8 


0.0375(7) 


0.0409(8) 


0.0409(9) 


12 


0.0357(17) 


0.0382(18) 


0.039(2) 


16 


0.0375(12) 


0.0395(12) 


0.0395(13) 


24 


0.038(5) 


0.038(5) 


0.035(6) 
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with the L\ = 8 curves. We obtain 



C = -2.4085(3), u = 0.86(4), x 2 /d-o.f. = 12.2/14, (33) 

that are fully compatible with a very slight increase of the statistical errors. This 
gives us confidence in the smallness of the higher order corrections-to-scaling. 

A further check can be done studying the crossing between lattice sizes with 
s = 2. The functional form of the fit is different now, as s in eq. (|3fj| ) is fixed. 
We obtain, discarding the L = 6 data, 

C = -2.4087(4), lu* = 0.83(4), x 2 /d-o.f. = 8.7/10. (34) 
where oj* is defined as u> + ^ — (0.783) -1 . 

We will use in the next subsection these results to study the infinite vol- 
ume extrapolation of the critical exponents. On the other hand, with our /3 C 
estimation, we can use again eq. (Gq) to obtain another estimation of critical 



exponents. This shall be done in subsection 5.2 



5.1 The infinite volume extrapolation of the exponents. 

The FSS analysis has the nice feature of using the finite-size effects to measure 
quantities as important as critical exponents. However, this does not guarantee 
the absence of systematic errors coming from finite-size effects, that appear as 
corrections-to-scaling terms. 

Let x be the critical exponent under study, and let us suppose that we need 
to take care of just the first term on the corrections-to-scaling series, then we 
have from eq. (p9|), 

- = - + aL-" + . . . , (35) 

V (2L,L) V 

a being a constant. 

We have carried out a fit, to the above functional form, for the numbers 



obtained in subsections 4.3.1 and 4.3.2. We exclude the L — 6 data, and take 



for the universal exponent to the value obtained previously (lu — 0.85(4)). 



Table 4: r\ exponent obtained from the susceptibility, magnetization and max- 
imum eigenvalue of the magnetization in (L, 2L) lattices. In the first row we 
indicate the quantity used as the variable. 



L 










kL 




X 


M 




X 


M 


\ max 


6 


1.442(2) 


1.447(2) 


1.442(3) 


1.332(6) 


1.332(6) 


1.332(5) 


8 


1.413(2) 


1.416(2) 


1.4135(18) 


1.332(4) 


1.332(4) 


1.334(4) 


12 


1.391(3) 


1.393(3) 


1.396(4) 


1.320(13) 


1.321(13) 


1.327(12) 


16 


1.381(3) 


1.383(3) 


1.384(3) 


1.339(5) 


1.340(5) 


1.343(5) 


24 


1.362(8) 


1.365(9) 


1.362(10) 


1.334(18) 


1.334(18) 


1.332(17) 
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Figure 5: Crossing points of scaling functions for pairs (Li,L 2 ), as a func- 



tion of L. 



for Li = 6 (solid lines) and L\ 



8 (dashed line), and 



L 2 = 12,16,24,32,^ 



We have found two very different situations. Some estimators present a non 
measurable size dependence (first four columns of table |2|, the whole table [| 
and the last three columns in table [l|). For them, we find no dependence on 
lo (when it moves within its error bars) of the extrapolated exponent, that has 
a value almost equal to the most relevant number in the tables (the L\ = 16, 
L 2 = 32 pair), but the errors are doubled. 

We also have found estimators with a quite measurable size dependence, 
namely the last two columns of table ^| (exponent v) , and the first three columns 
of table H (exponent 77) . The extrapolated values are 

u b = 0.800(9), is 6 = 0.804(8), 
771 = 1-338(6), 772 = 1.339(6), 773 = 1.346(6), [M) 

where the subindices refer to columns in tables || and ||. We find here an in- 
dependence of the extrapolated exponents of one half of the quoted statistical 
error. While the extrapolation for the 77 exponent is fully consistent with the 
quoted numbers in the last three columns of table f|, the v we get looks too big. 
Although both estimations are within two standard deviations, higher order 
corrections-to-scaling are apparent by plotting the last columns of table ^| as a 
function of . 

As a final result for the exponents, we choose operators that give stable 
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0.5 1.0 1.5 2.0 2.5 3.0 

CO 

Figure 6: Critical coupling extrapolated from the intersection with L\ = 6 for 
k, k Si £/L and £, S /L as a function of the lo exponent (upper part). Fits to a 
single (3 C value for L% = 6 (solid line), L\ = 8 (dashed line), and L2/L1 — 2 
(dotted line) as a function of u> (middle). In the lower side we show the values 
of x 2 f° r the latter fits. 

values with the lattice size, taking the measure for the (16, 32) pair; the number 
in the second parenthesis shows the error increase due to the extrapolation that 
can be taken as an estimation of the error involved in this procedure: 

v = 0.783(5)(6), 

r? s = 0.0380(12)(14), (ft = 0.406(3)(3), 7s = 1.536(10)(12) ), 
r) = 1.339(5)(5) (/3 = 0.916(6)(7), 7 = 0.518(6)(7) ). 

(37) 

For comparison, we quote the MC results for v obtained for the three di- 
mensional ferromagnetic 0(3) |30) and 0(4) |Q models: 

KO(3)) = 0.704(6), 

^(0(4)) = 0.7479(90). [6S> 

Our results for the RP 2 model are completely incompatible to those of the 0(3) 
one. Regarding the 0(4) model, although the differences are small, the value 
for the v exponent seems also incompatible. 



20 



5.2 Critical exponents measuring at (5 C . 

Taking the values j3 c = —2.4086(4) and u> — 0.85(4) we will repeat the mea- 
surements of the critical exponents using the measures of the usual operators 
at a fixed coupling. As in the determination of (3 C we use all the MC data, it 
is difficult to perform a correct computation of the statistical errors taking into 
account all correlations. In addition, there are several sources of systematic 
errors, consequently the evaluation of the errors on the exponent is harder. 

The functional form we fit is AL X / U (l + bL~ UJ ), where u> is fixed. We present 
the results for the computed exponents in table |^. In the first column we show 
the minimum value of L that we include in the fit. We remark that the inclusion 
of the L~ u term is essential to get a correct fit. To estimate the errors, we have 
extrapolated at our measured (3 C and at one standard deviation, adding the 
difference to the statistical error in the extrapolation. The u> error contribution 
is negligible. 

For the r\ exponents, as the slope of M and x as a function of (3 is large 
at the critical point, the errors in the determination of (3 C affect strongly that 
of rj. This error is, in part, overestimated since we do not take into account 
the statistical correlation between the measures used to determine (3 C and the 
magnetization operators. This, albeit complex, could have been done. However, 
as the f3 c value has been obtained after a somewhat elaborated procedure, the 
influence of systematic errors is difficult to quantify. Consequently we consider 
more solid the determination obtained from the quotient method. 

From the hyperscaling relation a = 2—dv, using ([5?]), we have a/v — —0.44(2) 
and therefore, no divergence is to be seen for the specific heat. We could never- 
theless stick to the FSS prediction. In fact in fig. [7] we show the specific heat at 
the critical point as a function of L a ' v . Strong scaling corrections are apparent, 
but the overall linear behavior is well satisfied. 

6 Vacuum structure in the critical region. 

The antiferromagnetic broken phase is characterized by a L — > 00 nonzero value 
for the magnetizations M and M s . However, studying the different terms of the 
tensors we can obtain a lot of information about the structure of the broken 



Table 5: Critical exponents obtained by measuring at f3 c . L stands for the 
minimum value of the lattice size we have used for the fit. The contribution of 
the error in j3 c and to is included. 



L 




V 






V 






d(/d/3 


Xs 


M s 




M 


6 


0.783(9) 


0.761(9) 


0.030(7) 


0.030(7) 


1.324(12) 


1.323(13) 


8 


0.784(15) 


0.771(14) 


0.023(9) 


0.023(10) 


1.316(16) 


1.316(16) 


12 


0.79(3) 


0.79(3) 


0.011(16) 


0.010(17) 


1.292(27) 


1.291(27) 


16 


0.75(7) 


0.82(8) 


0.02(3) 


0.02(3) 


1.31(5) 


1.31(5) 
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phase. We shall consider the behavior of several operators constructed from the 
magnetization tensors. As the staggered magnetization behaves as t@ B {(3 S ~ 0.4) 
and the magnetization as 2r (/? ~ 0.9) the former is dominant for t small, t being 
the reduced temperature. 

The first question we address is the breakdown of the even-odd symmetry. 
Near to the critical point it is very difficult to reach an asymptotic (L — * oo) 
behavior that allows us to know if that symmetry is broken, but instead we can 
use FSS arguments to find the answer. 

Let us consider the tensors product trace 

trM s M= i(trM2-trMjj), (39) 

where M c = e ven^> an d analogously for the odd sublattice. If the 

even-odd symmetry were not broken, the M S M tensor should go to zero as 
1/V = L~ 3 , while the LHS of eq. (§|) is bounded by a L~^' + ^ u behavior. In 
practice we have studied the quantity A = ((tr M S M) 2 ). 



We have used again the technique described in subsection 4.3 to determine 



the critical exponent. We obtain, after anL-+co extrapolation the value 

— = -3.389(15), (40) 

which is compatible with —2{[3+f3 s )/v = —3.377(13), and discards the —6 value. 
We thus conclude that the even-odd symmetry breakdown does occur. 

The naive picture that comes from the previous conclusion is that one of 
the sublattices is plane aligned and the other fully aligned, and the 0(2) global 
symmetry is preserved. To check this we have studied the conmutator of the M s 
and M tensors. If there were a remaining 0(2) symmetry, both tensors could 
be simultaneously diagonalized and the conmutator should vanish. In a finite 
lattice, according to FSS it should behave as L~ d ' 2 ~P e >" . Our results show that 
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Figure 8: Histogram of trM^ — trM^ at j3 = —3.0 for several lattice sizes. 



the exponent for B — (tr[M s , M] [M s , M]t) has strong corrections-to-scaling, 
but after a L — > oo extrapolation we obtain 



— = -3.406(11). 
v 



(41) 



Due to these corrections the error bar is probably underestimated. We observe 
that is again very different from —4.04(2) and compatible with — 2(/3 + (3 s )/v. 
The consequences of these results are that the 0(2) symmetry must be broken, 
and that the alignment directions of both sublattices are not orthogonal. 



7 Numerical results at low temperature. 

We present in this section some numerical results obtained from simulations 
in the low temperature (p < f3 c ) region. The very first thing to establish is 
the existence of a low temperature phase, where both order parameters are not 
zero, and where the symmetry between the even and odd sublattices is broken, 
as discussed in the previous section. We can answer both questions computing 
an histogram of tr(MM s ). This is shown in fig. || for (3 = —3.0. We see a 
clearly nonzero thermodynamic limit for tr(MM s ). For (3 < —3, the double 
peak structure is even clearer. 

We have next studied the second-neighbor energy for the odd and even sub- 
lattices independently. In fig. ^| we plot the histograms of the sublattice energy. 
We recall that for (3 = — oo, E^ = 1 for one of the sublattices and Ei — 0.5 for 
the other. We observe a clear double peak structure at (3 = — 4 (upper side). At 
= — 10 (lower side) the double peak structure is already clear even for lattice 
sizes as small as L = 4 and for L > 8 the peaks are completely split. The dotted 



line corresponds to the flat model (subsection 2.5), this means that for this (3 



value the results are almost asymptotic ((3 — > — oo) for the left peak. 
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It is very interesting to observe the evolution of the second-neighbor energy 
for the less aligned sublattice. At —f3 very large, the field variables are almost 
in a S 1 subset of S 2 , so, we can look at it as an RP 1 model in a face centered 
cubic (fee) lattice. This model (that is equivalent to an 0(2) sigma model) has 
a second order phase transition. For a rough comparison we have measured the 
energy (cos 2 (ipi — tpj)) for the fee RP 1 model at the critical point, obtaining 
E 2 = 0.62(1). 

In the RP 2 model, the fluctuations produce an effective coupling between 
second-neighbors of the less aligned sublattice. For — (3 very large, we can use the 
flat model defined by fl2C|). We obtain E 2 ~ 0.53 for the less aligned sublattice 
(see fig. ^). This value is very small, so the effective coupling is too weak to 
produce an 0(2) ordering. 

For increasing values of f3, the fluctuations off the plane reduce the E 2 value 
although (cos 2 ((fii — cpj)) increases. This effect can be taken into account in a 
crude approximation as follows. Let us write the spin variable as in eq. (|l9[). 
The second-neighbor energy takes the form 



((vi-vtf) = ((l-«) 2 )(l-( V |) 2 )cos 2 (^-^))+(«^ 



) 2 > 

(42) 
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As a first approximation, we can suppose that the z components of the neighbors 
are uncorrelated, and similarly for the z components with the angles. Thus we 
can factorize the mean values obtaining 

ig^McoB'te -¥>,-))« ^\ X h , (43) 

\L ~ ^min ) 

where A m i n = {(vf ) 2 ) is the minimal eigenvalue of the magnetization tensor of 
the less aligned sublattice. 

The numerical results obtained for several (3 values are 



= 


-10, 




= 0.53(1), 


(3 = 


-4, 


£ Rpl 


= 0.61(1), 


P = 


-3, 


£ Rpl 


= 0.65(1), 


(3 = 


-2.41, 


£ Rpl 


= 0.798(2) 



(44) 



We observe a fast increase of that quantity when approaching to the critical 
point. In particular, the critical value for the 0(2) model is reached at (3 = —4.0. 
This argument, although qualitative, supports the idea that the breakdown of 
the 0(2) symmetry observed at the critical point can be explained with the 
induced interaction of the more ordered sublattice on the less ordered one. 

We expect that the breakdown of the 0(2) symmetry occurs at (3 < (3 C since 
E RP becomes very large at f3 c . However, we have been unable to obtain a clear 
thermodynamic limit at (3 < (3 C for an 0(2) broken phase. Let us state that for 
an 0(2) symmetric phase, we expect to find a diminishing difference of the two 
smaller eigenvalues of the global magnetization as L~i . In fig. we plot this 
difference as a function of for two j3 values. At (3 = —4, the analysis for 

L < 32 lattices shows an 0(2) symmetric thermodynamic limit. At (3 = —3, the 
same analysis seems to indicate an 0(2) non-symmetric behavior for L < 32, 
but the L — 48 data points to a restoration of the symmetry. Finding a clear 
non symmetric limit would imply to simulate nearer to [3 C in much larger lattices 
than those we have used in this work. 

There is an alternative procedure that we think easier and more conclusive 
that consists in the introduction of a second-neighbor coupling which explicitly 
breaks the 0(2) symmetry 

5 2 = /3 2 Y, («i-«j) 3 - ( 45 ) 

<i,j> 

Let us discuss some regions of the (/3, /3 2 ) plane. For (3 = — oo, the line /3 2 
corresponds to a fee RP 1 (= 0(2)) model that presents a second-order phase 
transition at /3 2 = 0.415(1). At —j3 large but finite, we expect a second order 
line, in the same universality class, almost parallel to the (3 axis. As the fluctu- 
ations of the more aligned sublattice help in the 0(2) ordering, this line should 
get closer to the (3 axis and presumably crosses it before (3 C . We project to study 
this enlarged parameter space shortly. 
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Figure 10: Difference of the two smaller eigenvalues of M as a function of V * 
for L = 16, 32 at (3 = -4.0 (dashed line), and for L = 12, 16, 32, 48 at (3 = -3.0 
(solid line). 

8 Conclusions and outlook. 

We have found very interesting properties in the AF sector of the RP 2 spin 
model in three dimensions. It presents two magnetization operators with dif- 
ferent critical properties, what is related with an odd property of the vacuum 
in the critical region: a complete breakdown of the global 0(3) symmetry of 
the action. We have accurately measured the critical exponents, using a modest 
amount of computer time, by means of a FSS method that performs very good 
even in the presence of important corrections-to-scaling. The critical exponents 
obtained are clearly different from those of the ferromagnetic 0(3) model, and 
hardly compatible with the ones reported for 0(4). 

We project to extend these results to several related models, with general- 
izations such as, for instance, the inclusion of second-neighbor interactions or of 
vectorial interactions. Work is already in progress fl3l|] on the four dimensional 
version of this model. We are particularly interested in the relevance, in d = 4, 
of the properties found here. In addition, it can enlighten the mechanism by 
which the AF transition seems to belong to a new Universality Class in d = 3. 
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A Details of the Mean-Field computation. 



To compute the function Q p denned in ( jig ) we start from the Hamiltonians 
written in terms of the tensorial products r, = Vi ® Vi, 

H = tenTj, (46) 

<ij> 

Hp = - trnQe - trT *Q°> (47) 

i even i odd 

where, given the antiferromagnetic character of our model, it is rather natural 
to use a (tensorial) mean field for the odd sublattice, and an independent one 
for the even sites. In this approximation the partition function factorizes as 

Z p = z v J 2 z v J 2 , z c , = / dvexp(itrrQ e , ), (48) 
Js 2 1 

and the mean value of r can be written as 

(t . o ) p = M e , o (Q 0O ) = (z^r 1 / dv t exp( i tr t Q eo ) . (49) 

JS 2 1 

We thus obtain the quantity to minimize, G P , as a function of Q c , Q a , and 
the temperature T, 

^-G p = -Tlog(z c z ) -2dtrM M + trM Q + trM Q . (50) 

It is easy to show that Q p is invariant if we add to Q (Q ) an arbitrary 
multiple of the identity. Thus, we can shift to zero, for each field, one of its 
eigenvalues. 

Minimizing general tensorial fields would be very complex, so we shall re- 
strict ourselves to the case [Q c , Q ] = 0, where they can be simultaneously 
diagonalized. From now on, we shall work in a basis where both fields are di- 
agonal. Let us remark that an 0(2)-symmetric vacuum or a non-symmetric 
one, can be considered in this class. We shall write Q e , (Q ) as the diagonal 
matrix with eigenvalues (qi, 52, 0) ((^3, 54, 0)). Consequently, M Ci0 are trace-one 
diagonal matrices, with nonnegative eigenvalues. 

The values of the qi parameters that minimize Q p satisfy the equations 

= tr ({2dM c - Qo)^J , * = 1,2, (51) 

and analogously for i = 3, 4 interchanging odd by even. 

It is easy to check that Q c o = is a solution of ( [5l| ) for all values of T. 
Nevertheless, this could correspond to a maximum as well as to a minimum of 
Q p . In fact, one would expect this to be the minimum for large T, which would 
turn into a maximum for T small enough. 
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Eqs. ( |5l| ) are a set of quite complex nonlinear coupled equations. To go 
further, let us assume that the system undergoes a continuous transition at a 
critical point T c . We can thus linearize them as = 1 Aji(T)qi, in order to 
find a minimum close to Q . o = 0. To have a non trivial solution, we require 
det A{T) = 0, which happens at 

*-r-a- (52) 

This nontrivial minima should be in the directions in g-space pointed by the 
vectors belonging to the kernel of the matrix A, that is spanned by the vectors 
(1,0,-1,0) and (0,1,0 — 1). We shall then restrict ourselves to the subspace 
<7i = — ^3, q2 = — <74. Notice that there are three directions in this plane, where 
the ordered phase keeps an 0(2) symmetry: (1,1,-1,-1), (1,0,-1,0), and 
(0,1,0,-1). Expanding eq. ( |49| ) in powers of it is straightforward to see 
that the linear terms are opposite for M and M Q in this subspace, so, at first 
nonzero order 

M s = i(Mo-Mo) oc q h 

M = I(M C + M )-±1 oc qf, (53) 
E-\ = trM G M -i oc qf. 

To calculate the dependence of qi on T we need to consider higher order 
terms in Q p . We obtain, using polar coordinates in the (<?i, 52) plane 

Q p = a + b(T - T c )s 2 + c(T)s 4 + 0(s 6 ), (54) 

where s = r\/2 — sin 29, and 6, c(T) > 0. At leading order we also find tr M 2 oc s 2 
For T > T c , s = is a minimum of Q p , while when T < T c , s = is a maxi- 
mum and the real minimum is at s m in oc (T c — T) 1 / 2 . Therefore, the mean field 
prediction for the magnetic and specific heat critical exponents is 

f3 s = 1/2, = 1, a = 0. (55) 

Up to this order, the symmetry of the ordered phase cannot be determined 
since there are minima for Q p along all directions in the (qi, 52) plane. The sixth 
order term does break this degeneracy but favors an 0(2) symmetric ordered 
phase. This is in contradiction with our numerical data just after the transition. 



B Classical continuum limit. 

We shall derive the naive classical continuum limit for the RP 2 antiferromag- 
netic model. As we are interested in an 0(2) broken phase, for consistency we 
should add to the action ([!]) a ferromagnetic second-neighbor term that, as we 
shall see, does not change the continuum functional form. 

Let us first state that any integral over the S 2 sphere can be written as 
an integral over the SO (3) group. If u is an arbitrary element of S 2 , for any 
function /, 

/ dvf(v) = I dRf(Ru), (56) 

Js 2 JSO(3) 



28 



where dv is the invariant measure over the sphere and dR the Haar measure 
over the SO (3) group, both normalized. This result is easily generalized to a 
multiple integral. In particular the partition function can be expressed as 

Z = J {\{dv}j c -*(»i.-. «~) =jf ^ri di? ^ c- S{RlUu — RnUn \ (57) 

The vectors {ui, . . . ,un} should be chosen for Ri to be a smooth function of the 
spatial position in the {3 — > — oo limit, with the second-neighbor ferromagnetic 
coupling large enough to ensure the 0(2) breakdown. For a selection of the type 

u •= I Uc ' * even ' (58) 
* [ u Q , i odd, v ' 

the action can be written, in terms of the R variables, as 

s =-jj£ E E k^) • (iw«o)f 

i even q— ±1 /Lt=l,2,3 

+ y H X! X! [(-Ri«e) • {Rt+ q ^+puU c )f + (even <-> odd) . ^ 

i even q,p— ±1 *t,i>=l,2,3 

We need to take it c • u = in order R to be smooth. Thus we can substitute 
R by its Taylor expansion, and the action becomes at leading order 



S w constant H — / d 3 x [(R(x)u c ) ■ (d fi R(x)u )] 
aJ R3 y 

2(3' f 

/ d 3 a; ((dfj,R(x)u c ) ■ (d fl R(x)u c ) + (even <-> odd)) . 

a JR. 3 „ 



(60) 



a being the lattice spacing. 

If we choose, for example, the values u c = (1,0,0), u Q = (0,1,0), we can 
write [(Ru e ) ■ (d^Ruo)] 2 = [(R T d^R)^] , and given that R T d ll R is antisym- 
metric one can readily obtain, after some algebra, 

S= [ d 3 xy2[P(R T d^Rf], (61) 
where P is the diagonal matrix + + -£). 
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